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ABSTRACT 

We calculate the 7-ray albedo flux from cosmic-ray (CR) interactions with the solid rock and ice 
in Main Belt asteroids (MBAs), Jovian and Neptunian Trojan asteroids, and Kuiper Belt objects 
(KBOs) using the Moon as a template. We show that the 7-ray albedo for the Main Belt, Trojans, 
and Kuiper Belt strongly depends on the small-body size distribution of each system. Based on an 
analysis of the Energetic Gamma Ray Experiment Telescope (EGRET) data we infer that the diffuse 
emission from the MBAs, Trojans, and KBOs has an integrated flux of less than ^6 x 10~ 6 cm' 2 
s _1 (100-500 MeV), which corresponds to ~12 times the Lunar albedo, and may be detectable by 
the forthcoming Gamma Ray Large Area Space Telescope (GLAST). If detected by GLAST, it can 
provide unique direct information about the number of small bodies in each system that is difficult to 
assess by any other method. Additionally, the KBO albedo flux can be used to probe the spectrum of 
CR nuclei at close-to-intcrstellar conditions. The orbits of MBAs, Trojans, and KBOs are distributed 
near the ecliptic, which passes through the Galactic center and high Galactic latitudes. Therefore, 
the asteroid 7-ray albedo has to be taken into account when analyzing weak 7-ray sources close to the 
ecliptic, especially near the Galactic center and for signals at high Galactic latitudes, such as the ex- 
tragalactic 7-ray emission. The asteroid albedo spectrum also exhibits a 511 keV line due to secondary 
positrons annihilating in the rock. This may be an important and previously unrecognized celestial 
foreground for the INTErnational Gamma-Ray Astrophysics Laboratory (INTEGRAL) observations 
of the Galactic 511 keV line emission including the direction of the Galactic center. 
Subject headings: elementary particles — Kuiper Belt — minor planets, asteroids — Galaxy: bulge 
— cosmic rays — gamma-rays: theory 



1. INTRODUCTION 

The populations of small solar system bodies (SSSB) in 
the asteroid belt between Mars and Jupiter, Jovian and 
Neptunian Trojans, and in the Kuiper Belt beyond Nep- 
tune's orbit (often called also trans-Neptunian objects - 
TNOs) remain the least explored members of the solar 
system. A majority of the MBAs and KBOs have their 
orbits distributed near the ecli ptic with a FWHM of the 
order of 10° in ecliptic latitude ^Binzel et al.l ll999: Brown 
2001). The spatial and size distributions of these objects 
provides important information about the dynamical 
evolution of the solar system. Extending our knowledge 
of the size distribution of these objects below current sub- 
kilometer size limits of optical (e.g., Ilvezic et~aLl 120011 : 
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Wieeert et al.ll2007f) and in frared (e.g-. lTedesco fc DesertJ 
2002t lYoshida et al.ll2003f) measurements would provide 
additional information on the accretion/collision and de- 
pletion processes that formed the populations of SSSBs 2 . 
In this paper we show that the CR-induced 7-ray albedo 
of these systems may be bright enough to be detected 
with a 7-ray telescope such as GLAST and/or INTE- 
GRAL and/or Soft Gamma-ray Detector ( SGD) aboard 
the NeXT satellite (jTakahashi et all 12009) (see our es- 
timates below), and can allow us to probe the size dis- 
tribution of SSSBs down to a few metres. Additionally, 
the 7-ray emission of these systems may comprise a "dif- 
fuse" 7-ray foreground that should be taken into account 
when evaluating the flux and spectra of 7-ray sources 

2 We note that Babich et al. (2007) have suggested a method to 
place constraints upon the mass, distance, and size distribution of 
TNOs using spectral distortions of the CMB. 
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ne ar the ecliptic. Our p reliminary results are presented 
in IMoskalenko et al.l (|2008l ). 

The Galactic center is a region crowded with 7-ray 
sources and is one of the preferred places to look for 7-ray 
signatures of dark matter (DM ). An extensive li t eratur e 



on the subject exi st s , 
IZahariias fc Hoopeil 



Bergstrom et alJ (fl 998). 



Finkbeiner fc Weinerl 
^ l200l ~ 



iBaltz et all (|2008D ; also 
The ecliptic crosses the 
Galactic center almost 



( 20071 ). IHooper et alJ (pOO 
references in these papers. 
Galactic equator near the 
perpendicularly with inclination ^86.5°, and underes- 
timation of the SSSB albedo foreground may lead to 
errors in the analysis of weak or extended sources in this 
region. 

The Galactic center region also harbors the enigmatic 
source of the 511 keV positron annihilation line observed 
by the Orien ted Scintillation Sp ectrometer Experiment 
(OSSE) (e.g-lPureell et al.lll997f) and INTEGRAL (e.g ., 
iKnddlseder et alJ 120051 : IWeidenspointner et all 120061 ). 
The distribution of the annihilation line does not match 
the distribution of any positron source candidate, e.g., 
pulsars, supernova remnants, binaries, radioactive iso- 
topes, such as 26 Al, etc. A number of excellent dis- 
cussions on the origin of this emission are available 
in the literature, ranging from positron focusing by 
the r egular Galactic magnetic field to DM annihila- 
tion (Finkb cincr & Wcincr 2007; Guc ssoum et al.l [2005: 
IHooper et al.ll2008t Uean et al.l l2006: Pra ntzosll2006l and 
references therein). Our calculations (detailed below) in- 
dicate the SSSB CR-induced albedo spectrum should ex- 
hibit a 511 keV line due to secondary positrons annihi- 
lating in the rock. Since the target material (rock, ice) is 
solid, the line has to be very narrow. This emission pro- 
duces a previously unrecognized celestial foreground to 
the 511 keV flux including the direction of the Galactic 
center. 

At higher energies, above ~30 MeV, regions at high 
Galactic latitudes are conventionally used to derive the 
level of the extragalactic 7-ray emission by compar- 
ing a model of the diffuse Galactic emission to the 
point-source-subtrac ted skymaps and extrapolating to 
zero m odel flux (e.g. JSreekumar et al. 1998; S trong et al.l 
2004). The remainder is assumed to represent the 
level of the isotropic, presumably extragalactic emission. 
However, recent studies have predicted another impor- 
tant foreground component with a broad distribution 
on the sky originating from the inverse Compton scat- 
tering of solar photons by CR electrons in the helio - 
sphere (jMoskalenko et al.ll2006l : IOrlando fc Stronsll2007D . 
which has to be included in the analysis of the diffuse 
emission. A reanalysis of the EGRET data revealed 
this broad componen t, in agreement with the predictions 
(|Orlando et al.ll200l . Since the ecliptic passes through 
high Galactic latitudes, the SSSB albedo flux also may 
need to be taken into account when analysing the weak 
extragalactic component. 

2. SMALL SOLAR SYSTEM BODIES 

The asteroid mass and size distributions are thought to 
be governed by collisional evolution and accretion. Col- 
lisions between asteroids give rise to a cascade of frag- 
ments, shifting mass toward smaller sizes, while slow 
accretion leads to the growth of the latter. The first 
comprehensive analytical description of such a collisional 



cascade is given by DohnanyJ (|1969f ) . Under the assump- 
tions of scaling of the collisional response parameters and 
an upper cutoff in mass, the relaxed size and mass dis- 
tributions approach power-laws: 



dN = am k dm 
dN = br~ n dr, 



(1) 
(2) 



where m is the asteroid mass, r is the asteroid radius, and 
a, b, k, n are constants. These equilibrium distributions 
extend over all size and mass ranges of the population 
except near its high-mass end. The constants in eqs. ([T]), 
@ are not independent. If all asteroids have the same 
density p, one obtains n — 3 k — 2 and b = 3a(AiTp/3) 1 ~ k 
(see eqs. [3], [4]). For a pure iDohnanvil cascade k = 11/6 
and n = 3.5. 

However, collisional response parameters are not size- 
independent, e.g., the energy per unit target mass deliv- 
ered by the projectile required for catastrophic disrup- 
tion of the target (the so-called critical specific energy) 
depends on the radius of the body, and the single power- 
laws (eqs. [Tj, [2]) break. Even though the sizes of aster- 
oids generally can not be directly observed (except by a 
small number of asteroids studied by spacecraft flybys, by 
stellar occultation, or those well observed by radar) and 
are instead estimated using apparent magnitude, optical 
and infrared albedos, and distances, the information col- 
lected on a large sample of MB As 3 seems to confirm that 
the real distribution departs from a single power law, at 
least for objects larger than a few kilometers. Smaller 
sizes are very difficult to detect, and one has bear in 
mind the observational bias of the incompleteness of the 
small (dim) ast eroid sample. Though de- biasing can be 
attempted (e.g.. lJedicke fc Metca lfe 1998), a large ambi- 
guity still remains. 

Figure[T]shows the MBA size distributions as published 
in the literature and those used in this paper. For the 
MB As larger than diameter V (km) , iBinzel et alJ (|1999l ) 
give N(> v) = 1.9 x 10 6 £>~ 2 - 52 (the authors do not give 
the range of sizes, so we adopted a cut at V ~ 0.5 km). 
ITedesco fc Desert! (|2002ft give \ogN(> v) = (5.9324 ± 
0.0016) - (1.5021 ± 0.045) log£> for 0.2 km < V < 2 km 
based on Infrared Space Observatory (ISO) observations. 
Using a sample of m ore than 6 x 10 4 MB As to a limiting 
magnitude of 7~21, Uedicke fc Metcalfe! (|1998h found a 
change in the slope of the cumulative distribution from - 
2.25 for 1 km < V <10 km to -4.00 for 10 km < V < few 
10s of km. Based on observations of ^13000 MBAs by 
the Sloan Digital Sky Survey (SDSS). IIvezic et al.l (pOOlh 
found that the cumulative size distribution resembles a 
broken power-law, oc I? -2 - 3 for 0.4 km < r < 5 km, 
and oc 2?~ 4 for 5 km < r < 40 km, a nd is indepen- 
dent o f the heliospheric distance. Finally, ITedesco et al.l 
(2005) gives a fit to data between 1 km < T> < 100 km, 
log AT (> V) = 6.275 ± 0.013 - (3.214 ± 0.056) log V + 
(0.974 ± 0.066) log 2 V - (0.182 ± 0.022) log 3 V, but ex- 
trapolation to smaller sizes is invalid. The size distri- 
bution below ^1 km is essentially unexplored territory. 
One piece of evidence comes from the size distribution of 
ejecta blocks on 433 Eros. Based on the block d istribu- 
tion over a size range 0.1 - 150 m, I Cheng] (|2004l ) argued 

3 The Minor Planet Center supports a database for all observed 
SSSBs: |http://www.c fa,harvar dTedu/iau/mpc.htmlj 
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Fig. 1. — Cumulative size distribution 2V(> r) of KBOs (upper 
set of lines ) and MBAs ( l ower s et of lines) . Line c oding: thick dash- 
dot li ne — [Binzcl ct al. ( 19 93), thick dash es — [Tcdcsco fc Desert] 
ll2002h , thick solid - ITedesco et al.l f2005h . thick dots - parame- 
terization proposed by Cheng (2004). Our parameterizations are 
shown by thin lines (solid, dotted) where the numbers show the 
cumulative power-law index (n — 1) of a particular distribution. 
Thin solid lines are our adapted distributions: index 2.0 (n = 3.0) 
for MBAs, and 2.5 (n = 3.5) for KBOs. Thin dotted lines show 
the range discussed in the paper. See text for details. 



that these data support a cumulative index 2.5 extrapo- 
lation down to sizes ~lm. Our distribution with a single 
cumulative index (n — 1) = 2 (thin solid line), detailed 
in the next Section, seems to match the global size dis- 
tribution determined from various types of observations 
in the wide range of radii 10 2 — 10 7 cm. We will use this 
distribution in our estimates of the MBA albedo, below. 

The dynamical estimate of the total mass of the aster- 
oid b elt is about (3.6 ± 0.4) x 10 24 g (jKrasinskv et all 
2002) or close to 5% of the mass of the Moon. The total 
mass is dominated by large bodies, while the 7-ray albedo 
is dominated by very small bodies. The largest MBA, 1 
Ceres, comprises about 30% of the total mass of the as- 
teroid belt alone. However, it does not provide a restric- 
tion on the size (and mass) distribution of small bodies. 
Current estimates indicate the total numb er of MBAs 
above 1 km in diameter is (1. 2 — 1.9) x 10 6 (iBinzel et al.l 
Il999t ITedesco fc Desert! l2002t ITedesco et al.ll2005ft . Our 
adopted distribution with n = 3 gives a number near the 
upper end of this range, 1.92 x 10 6 , while also putting 
the total number of MBAs with r > 1 m at ~ 5 x 10 11 
(Figure [J) . To get an idea of how the MBA albedo flux 



depends on the extrapolation to small radii T> < 1 km, 
we also consider broken power-law distributions with in- 
dices 2.5 and 3.5 below 1 km in diameter, retaining an 
index 3 for sizes larger than 1 km. 

The densities of most M BAs lie in the range 1.0 - 3.5 g 
cm -3 ( Bi nzel et al.lll999l ) while the densities of particular 
asteroid classes can vary broadly, 1.23 - 1.40 g cm' 3 for 
carbonaceous, 2.65 - 2.75 g cm~ 3 for silicate, and 4.75 - 
5.82 g cm -3 for metallic bodies (jKrasinskv et al.l l2002). 
We adopt an average density p = 2 g cm~ 3 . 

Most MBAs have a semimajor axis between 2.1 and 
3.3 AU with a low eccenticity orbit. In our estimates we 
assume an average circular orbit with radius £ ~ 2.7 AU. 

The Jovian Trojan populations of asteroids are collec- 
tions of bodies in the same orbit as Jupiter (semimajor 
axis £ ~ 5.2 AU) located at the L4 and L5 Lagrange 
points of the Jupiter-Sun system. The Trojans are thus 
concentrated in two regions rather than distributed over 
the entire ecliptic as for the MBAs. The total mass of the 
Jovian Trojans is estimated to be ~ 10 -4 M© where Mq 
is the mass of the Earth with a differential power-law in- 
dex n ~ 3 in the size range 2 km to 20 km (jJewitt et al.l 
l2000UYoshida fc Nakaniurall2005| ). similar to MBAs, giv- 
ing a number of ob jects > 1 km in diameter ~ 1.3 x 10 6 
(jJewitt et al. 2000). The combined mass of these objects 
is approximately the same as for the MBAs. The number 
of objects > 1 km in diameter and the power-law index 
n ~ 3 makes their size distribution very similar to that 
of MBAs. 

The mass density of SSSBs in this group varies sig- 
nificantly: estimates for the binary Trojan 617 Patroclus 
are less than water ice p = 0.8to' 2 g cm -3 (|Marchis et al.l 
[20061) . as are those for other Trojan binaries p ~ 0.6 — 0.8 
g cm -3 (|Mann et al.ll2007t ), while 624 Hektor is some- 
what denser p = 2.48^o'ggQ g cm -3 (jLacerda fc Jewittl 
I2007I) . In our calculations we adopt an average density 
p = 1 g cm -3 as a compromise between these bounds. 

We also consider ic y bodies and comets in the Kuiper 
Belt (for a review see ILuu fc Jew itt 2002|) and the con- 
joining innermost part of the Oort Cloud 4 , but call them 
all KBOs for simplicity. The KBOs are not uniformly 
distributed, with at least three dynamically distinct pop- 
ulations identified: the Classical Disk, the Scattered Disk 
with large eccentricities and inclinations, and "Pluti- 
nos" around the 3:2 mean motion resonance with Nep- 
tune at 39.4 AU. Kui per Belt Objects are distributed 
between 30 - 100 AU (|Backman et al.lll995l and refer- 
ences therein) with surface number densit y a(£) — A£~ a 
(jBackman et al.lll995UJewitt fc Luulll995| ). where A is a 
constant determined in eq. (TTTj) , and a — 2. The total 
mass is estimated to be in the range ^0.01 -0.3 M ffi , while 
the m ost often used value is ~0.1 -M© (|Luu fc JewittJ 
l2002h . The density of small icy bodies and com ets is 
~0.5 g cm" 3 (lAsphaue fc Benzlll99l : lSoleml[l99l . 

The KBO size distribution is much more difficult to 
determine because of their dimness. It is widely believed 
that the TNOs are dynamically related to the Centaurs 

4 The Oort Cloud of comets (e.g., Stern 2003) is thought to oc- 
cupy a vast space between 50 and 50000 AU from the Sun and also 
contributes to the celestial 7-ray foreground. However, its exact 
mass and distribution are poorly constrained. We are planning to 
investigate limits on the albedo of the Oort Cloud in a forthcoming 
paper. 
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(planetesimals distributed between Jupiter and Neptune 
that are in crossing orbits of the giant gas planets), and to 
the Jupiter-family group of ecliptic comets that may be 
objects that were knocked inwards from the Kuiper belt. 
The KBO size distribution is determined by very indirect 
methods such as measuring the sizes of the nuclei of the 
ecliptic comets (and making assumptions on how they 
evolve during their repeated passages through the inner 
solar system), and Centaurs, as well as impact craters 
on the Galilean satellites of Jupiter. The estimates of 
the size distribution for the cometary nuclei range from 
n = 2.6 to 3.7 in the range r ~ 1 — 10 km, for KBOs - 3.7- 
4.45 (r > 20 km), and for Centaurs - 3. 7-4.0 (an appro- 
priat e discussion c an be found, e . g., in iBernstein et al.1 
12004 iTothl 120051 . iTancredi et all 120061 and references 
therein). If the index is > 4.0, the mass of the total 
population is dominated by the smallest bodies. How- 
ever, there are some reasons to believe that the size dis- 
tribution begins to flatten wel l above 1 km in size. Col- 
lisional evolution simulations (jKenvon fe Bromlevll2004t ) 
show that the size distribution is a power law with index 
~4.5 for large bodies (r > 10 — 100 km) and — 3.5 - 4 
for small bodies (r < 0.1 — 1 km) for a wide range of 
bulk properties, initial masses, and orbital parameters. 
Adopting a conservative value of n = 3.5, we obtain the 
total number of comets (V > 1 km) at ~ 9 x 10 9 , whi ch 
is in agreement with other estimates (e.g.. ISternl 12003). 

There are als o large pop ulations of Centaurs NiV > 
1 km) - 10 8 (|Sheppardl l2"005T l between the orbits of 
Jupiter and Neptune, and Neptunian Trojans at the L4 
and L5 points of the Neptune-Sun system. The number 
of large Neptunian Trojans (V > 80 km) outnumbers 
the number of large Jovia n Trojans by a factor ~ 10 
dSheppard fc Truiilld[2006T ). Their power-law index may 
be close to that of the KBOs n ~ 3.5 making their 7-ray 
albedo essentially brighter than MBAs and Jovian Tro- 
jans at the same distance. While Centaurs are scattered 
between Jupiter and Neptune, the positions of Neptu- 
nian Trojans are well known so that the detection of a 
7-ray albedo signal may be simplified. 

3. CALCULATIONS 

We use the Lunar albedo spectrum as an approxima- 
tion of the SSSB albedo for two main reasons: (i) the 
Moon is a solid body in which the CR cascade in the 
rock develops similarly, and (ii) its proximity to the 
Earth allows it to be easily detectable by 7-ray tele- 
scopes. The spectru m of 7-rays from the Moon ha s been 
calculated recently (|Moskalenko fc Porter! l2007al fbl) us- 
ing the GEANT4 Monte Carlo framework to simulate 
the CR cascade development in a Lunar rock target (re- 
golith). It has been shown that the Lunar albedo spec- 
trum is very steep with an effective c utoff around 3-4 
GeV in agreement with observations (|Thompson et al] 
1997). The central part of the disk of the Moon has an 
even steeper spectrum with an effective cutoff at ^600 
MeV. The emission above ^10 MeV is thus dominated 
by the thin rim where CR particles interact close-to tan- 
gentially with the surface and the high-energy secondary 
7-rays shower out of the Moon toward the observer. In 
contrast to other astrophysical sources, the albedo spec- 
trum of the Moon is well understood, including its ab- 
solute normalization; this makes it a useful template for 
estimations of the CR-induced albedo of SSSBs without 



an atmosphere. Since the Moon functions as a standard 
(7-ray) candle, in our estimates we use the flux of the 
Moon as our standard and introduce the term "Lunar 
albedo flux units." 

If the SSSB size distribution dN/dr is known, it can 
be directly substituted into eq. ([6|) to estimate the 7-ray 
albedo flux. Below, we derive this albedo flux assuming 
that the size distribution for a single SSSB population 
is a simple power law where the normalization has to be 
obtained from the total mass of the system. 

Let the SSSB mass distribution have the form given by 
eq. ([T]) , which can be rewritten as a size distribution (cf. 
eq. 0) 



dN dN dm 



,dN 



dr dm dr ^ dm ' 



(3) 



where p is the average density of the SSSB target. As- 
suming that all SSSBs have a spherical shape, m = 
(4n/3)pr 3 , we get 



dN /4tt 



i-fc 



„2-3/c 



The normalization for a is obtained from 

1 dN 
m—dm = /Mmooii, 
dm 



(4) 



(5) 



where m and mi are the lower and upper SSSB masses, 
respectively, and /Mmooii is the total mass of the SSSB 
emitting population considered as a fraction / of the 
Moon's mass (/ = 0.05 for MBAs, / = 0.1M e /M M oon ~ 
8.13 for KBOs). The flux of 7-rays from such an ensem- 
ble of bodies with size distribution dN/ dr is then 



F = F Mo 



'Mc 



d J J rn dr Rmo 



-dr 



(G) 



-RMoon V d / V 3 / 



where 



ri 



1 



Un(ri/r ), 

fc-2 



k = 4/3 



47T 



a = /Mm 



....11 \ v / 'I 



2-k 

6-3fc 6-3fc ' r 2 



iln 1 {r 1 /r Q ), 



' 



k = 2 



Here -Fmooh is the Lunar rim albedo flux for the same 
incident spectrum of CR particles, i?Moon = 1-7382 x 
10 s cm is the Lunar radius, £>Moon — 0.0025 AU is the 
Earth-Moon distance, d is the distance (in AU) of the 
SSSB population from Earth, and ro and r\ are the sizes 
corresponding to the masses mo and mi . 

The factor r/Ry[ oon in eq. © comes from the fact 
that the albedo of SSSBs is dominated by the emission 
from the rim. The rim and the disk albedo fluxes of the 
Moon are about equal at low energies (Figure[2]), with the 
rim albedo flux considerably dominating above 10 MeV 
(|Moskalenko fc Porterir2007bl ). Since the rim albedo flux 
scales tx r, and the inner part of the disk oc r 2 , as the 
size of the emitting body decreases it is the rim which 
produces most of the albedo photons for SSSBs. 
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Assuming k ^ 4/3, 2, after some rearrangement we 
obtain 



3 

F{£>8) = TJ-Z^Moon^Moon 
Z7T 



PMo 



D 



Moon 



G(ri,r ;fc) 



1 



4-3fc 



1 



d(£,9)_ 
(r /n) 4 - 3fc 



l - (ro/n)^ 



G(ri,r ;k), 
(7) 
(8) 



where p M oon 
Moon, 



3.3 g cm 3 is the mean density of the 

iV2 



d(e,0) = cos9+(£ 2 



sm 



(9) 



£ (AU) is the radius of the orbit of the SSSB population 
(for MBAs, £ = 2.7 AU; for KBOs, see below), 9 is the 
angle between the line of sight (in the ecliptic) and the 
direction to the Sun, and we divided by 2tt to obtain the 
flux per radian. The total flux integrated over 9 is 



Ftat — 



F{£,6)d9. 



(10) 



For the case of the KBOs distributed between £ — 30 and 
100 AU, an additional integration over £ is required 



2tt 



100 



f£t = / d9 / F{£,9)a{£)£d£, (11) 

JO J30 

where a{£) = At~ 2 and A = -//ln0.3. These formulae 
provide fluxes integrated over ecliptic latitude. 

For the case of MBAs, our adopted differential size 
distribution is a broken power law with index ri\ = 3fci— 2 
for r > rt, and n% — 3^2 — 2 for r < r&: 



dN *tvr f 4n 

—j— — 6 J JWMoon I -g-p 



2-ki 



(12) 



-2-3ki 



r>r b 



6-3/ci _6-3fci 1 r 2-3fci 
r l ~ r 1 U 



n 



2~3k 2 



,r <r h 



where rt, = 0.5 x 10 5 cm, and we assume 2 < < 4. 
This can be inserted into eq. ^ to derive corresponding 
expressions for the flux from such a population of SSSBs. 

We can see that the observed albedo flux gives di- 
rect information on the integral / drr(dN/dr), eq. (J6]), 
which can be used to constrain the effective average ra- 
dius of the emitting bodies (r) and their total number 
in the system. Additionally, if the size distribution is 
a single power law, the observed albedo flux can pro- 
vide us with information about the power-law index. As 
can be seen from eq. (|8|). the function G(r±,ro; k) is a 
steep function of k. For k < 4/3, the expression in the 
last square brackets (eq. [S]) is ~1 since ri » r and 
G w rf 2 (2 - fc)/(4 - 3fc). For k > 2, eq. © becomes 

G ~ r o^ 2 (2 — fc)/(4 — 3fc). In the intermediate region 
4/3</c<2,G quickly increases with k. For the distri- 
bution of radii we consider the corresponding range for 
the size distribution index is n = 2.5 — 3.5 for MBAs 
and 3.0-4.0 for KBOs. This translates into an index, k, 
for the mass distribution in the range 1.5-2.0. In this 
range G{r\, r$; k) changes by 3 - 5 orders of magnitude 
depending mostly on the assumed value of rg. This al- 
lows a determination of k assuming the average density 



of the asteroid rocks is known. The function G also con- 
tains a dependence on r% , the radius of the largest body. 
For the MBAs, we use Ceres, r% = 4.565 x 10 7 cm, and 
for the KBOs we use 134340 Pluto, r x = 1.195 x 10 8 
cm, but the exact value of n does not change the size 
distribution significantly and does not affect our conclu- 
sions. Figure [1] shows the adopted size distributions of 
MBAs and KBOs which agree well with those given in 
the literature. 

The question of where most of the heliospheric mod- 
ulation occurs is important for the determination of the 
CR flux at an arbitrary distance from the Sun. The re- 
cent crossing of the heliospheric ter mination shock b y 
the Voyager 1 spacecraft at ~94 AU ([Stone et al.ll2005h . 
currently at ^104 AU, while Voyager 2 is still inside the 
termination shock, allows unique studies of the spectra 
of CR particles on both sides of the shock. Low-energy 
CR detectors on board the spacecraft indicate that the 
particle spectra are significantly different, supporting the 
conclusion that a considerable modulation of the CRs 
occurs near the termination shock. On the other hand, 
most of the albedo emission discussed in this paper is 
produced by CR particles with energies >1 GeV; their 
flux does not change significantly from local interstellar 
space down to ~40 AU, as indicated by current helio- 
spheric models (e.g.. lLangner et "aLll2006D . 

The Lunar albedo flux , Fmooth is calculated u sing th e 
procedure described in [Moskalcnko fc Porter! (|2007bD . 
To calculate the Lunar albedo at an arbitrary modu- 
lation level, we use the local interstellar (LIS) spectra 
of CR protons, helium, and positrons, as fitted to the 
numerical results of the GALPROP propagation model 
(jPtuskin et alj 2006, Table 1, re acceleration and plain 
diffusio n models) as described in iMoskalenko fc Porter! 
(|2007bl eq. (2) with parameters listed in Table l) 5 . The 
CR particle flux at an arbitrary phase of solar activity 
can t hen be estimated using the force-field approxima- 
tion eeson fc Axfordl[T96l : 



dJ p {E k ) _ dJ™(E k + &Z/A) E 2 



M 



J2 



dE k 



dE k 



{E + <S>Z/A) 2 - M 2 ' 



(13) 

where dJ£° jdE k is the LIS spectrum of the CR species, 
E k is the kinetic energy per nucleon, E is the total energy 
per nucleon, <!> is the modulation potential, Z and A are 
the nuclear charge and atomic number correspondingly, 
and M is the nucleon mass. The modulation potential 
$(£) at an arbitrary distance £ from th e Sun can be calcu- 
lated using the expressions derived in IMoskalenko et all 
|20p), their eqs. (7), (8). 

Figure [2] shows the Lunar albedo spectrum for differ- 
ent modulation potentials $ = 0, 500, 1500 MV. The no 
modulation case (<& = 0) corresponds to the upper limit 
of the KBO albedo, with moderate modulation ($ = 500 
MV) corresponding to the MBA albedo. The difference 
in brightness below ~1 GeV due to the incident CR flux 
only (no modulation vs. moderate modulation) is as large 
as a factor of ~2-3. Also shown are the components of 
the albedo spectrum (center, rim) for $ = 1500 MV. 

Figure [3] shows the albedo spectrum of the Moon as if 

5 The parameterizat ion constants for CR positrons, not given in 
Table 1 of IMoskalenk o fc Porter] l!2007bO . are: J = 44.8143, <*i = 
1, 6i = 0.594634, ci = -9.14888, a 2 = -605.291, b 2 = 1.53611, 
c 2 = -7.27809, a 3 = 1.18135, b 3 = 0.365787, c 3 = -3.51576. 
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Fig. 2. — Calculated 7 -ray albedo spectrum for CR nuclei inter- 
actions in the Moon rock Moskalenko & Porter (2007b) for selected 
modulation potentials. Line colouring: black, no modulation; red, 
$ = 500 MV; blue, $ = 1500 MV. Dashed and dotted lines show 
the albedo of the disk and the rim correspondingly for the case of 
<S> = 1500 MV. 



the Lunar surface was made of different materials: water 
ice (multiplied by a factor of 0.1), regolith, and iron (mul- 
tiplied by a factor of 10). The plot shows the albedo spec- 
tra for two limiting cases, no solar modulation ($ = 0) 
and solar maximum conditions at 1 AU ($ = 1500 MV). 
The low-energy parts of the spectra (<10 MeV) from 
different materials are considerably different and the nu- 
clear emission lines can be used to distinguish between 
the materials. The high-energy parts are essentially fea- 
tureless and have similar shape. The flux between iron 
and water ice changes by a factor of ^2 above 100 MeV 
with the latter producing the larger flux. Above ^100 
MeV the regolith albedo approaches the water ice albedo. 

The 511 keV line in Figures [2] and [3] is due to the an- 
nihilation of positrons produced by CR cascades in the 
solid target (iron, regolith, ice). In Figure [5J the albedo 
spectrum also includes the contribution by CR positrons 
in the Lunar rock target (see below). Since the rock is 
solid, secondary positrons quickly thcrmalize and pro- 
duce a narrow annihilation line. Its width is determined 
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Fig. 3. — Calculated 7-ray albedo spectrum of a Moon-sized body 
at the Lunar distance composed of moon rock (black), iron (xlO, 
blue), or water ice (xO.l, red). Line-styles: solid, no modulation; 
long-dashed, <t> = 1500 MV. 



by the energy bin size adopted in the simulation. 

Figure 0] shows the components of the albedo spec- 
trum (Figure HJ) below 20 MeV. The thick solid lines 
are the total albedo flux due to the CR proton, helium, 
and positron interactions with regolith for no modulation 
(upper, red) and modulation level 1500 MV (lower, blue). 
The thin solid lines show the albedo spectrum due to CR 
positron interactions with regolith for the same cases of 
no modulation (upper, red) and modulation level 1500 
MV (lower, blue). The dashed and dotted lines show the 
components of the CR positron induced 7-rays, from the 
center and the rim, correspondingly. 

4. ANALYSIS OF THE EGRET DATA 

The EGRET instrument on the Compton Gamma- Ray 
Observatory (1991-2000) surveyed the sky in the range 
>30 MeV and here we use the EGRET data together 
with the i nformation in the T hird EGRET Source Cata- 
log f3EG. lHartman et al.lll999l ) to set limits on the signal 
from the SSSBs. Challenges to detecting diffuse emission 
associated with the ecliptic plane in clude the brightnes s 
of the Galactic diffuse emission (e.g-. lHunter et all 1997). 
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Fig. 4. — Calculated 7-ray albedo spectrum showing components 
below 20 MeV for no modulation (red) and modulation level 1500 
MV (blue). Line-styles: long dash: positron induced 7-rays from 
center; short dash: CR positron induced 7-rays from rim; thin 
solid: total CR positron induced 7-rays; thick solid: total 7-ray 
emission from CR positrons and nucleons. 



the presence of bright point sources, the limited angular 
resolution and photon statistics of the EGRET data, and 
potential large-scale artifacts in the exposure maps owing 
to ageing of the spark chamber gas. 

We made maps of the EGRET data in ecliptic coordi- 
nates for Cycles 1-4 of the mission, during which most 
of the EGRET exposure was obtained. The event data, 
after standard cuts on zenith angle and inclination angle 
(30°), were binned on a photon- by-photon basis in eclip- 
tic coordinates. The exposure maps for each EGRET 
viewing period were transformed into ecliptic coordinates 
and added together and intensity maps were calculated 
from the photon and exposure maps. 

In order to limit contributions from Galactic diffuse 
emission to any possible enhancement of diffuse inten- 
sity at low ecliptic latitudes, the region \b\ < 20° for 
\l\ < 90° and \b\ < 10° for \l\ > 90° was masked out in 
the analysis. We also removed regions 12° in diameter 
around the position of each identified source in the 3EG 
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Fig. 5. — Profiles of 7-ray intensity with /3 derived from EGRET 
data as described in the text. The energy range is 100-500 MeV 
and the profiles have been averaged over all ecliptic longitudes, (a) 
Profile derived with no masking of Galactic diffuse emission or 7- 
ray point sources, (b) Profile with the Galactic plane (|6| < 10° 
for \l\ > 90° and \b\ < 20° for \l\ < 90°) excluded , (c) Profile with 
the identified 3EG sources (Hartman ct al. 1999) and the Galactic 
plane excluded, (d) Profile with the identi fied 3EG sources p lus the 
further blazar identifications proposed by Sowards-Emmerd] ( 12003 . 
2004) excluded. Overlaid on each profile is the best-fitting gaussian 
(12.5° FWHM, centered on j3 = 0) plus a constant, fit for the 
region \/3\ < 50°. This approximates the distribution of albedo 
7-ray emission expected for the the KBO. 



catalog (designation other than "u" in the 3EG catalog). 

Figure [5] presents the profile of 7-ray intensity in the 
100-500 MeV range over ecliptic longitude. This range 
was chosen as having the brightest expected albedo emis- 
sion in the energy range of EGRET. As described in the 



8 Moskalenko et al. 



Table 1. Diffuse intensity around the ecliptic (100-500 MeV) 



Set of cuts 


\P\ < 15 


o 








in Figure [5] 


Flux, cm 2 s 1 


Stat, error 


Fitted flux, cm 


" 2 S" 1 


Stat, error 


a 


1.006 X 10" 5 


5.5 x 10" 7 


9.16 X 10" 


6 


3.5 X 10" 7 


b 


7.95 x 10-° 


5.8 x 10~ 7 


5.95 x 10- 


6 


3.7 x 10~ 7 


c 


3.59 x 10~ 6 


6.7 x 10~ 7 


3.53 x 10- 


6 


4.4 x 10~ 7 


d 


1.1 x 10" 7 


7.4 x 10~ 7 


1.52 x 10" 


6 


5.1 x 10" 7 



caption, a sequence of profiles is shown for different com- 
binations of the masks described above. In the last pro- 
file (Fu2ure^)j_thj3E^^urces subsequently identified 
bv lSowards-Emmerdl (|2003L l2004h as likely to be blazars 
were included with the sources identified in the 3EG cat- 
alog in defining the mask. We did not mask out uniden- 
tified point sources because of the possibility that some 
of them may have represented detections of the 7-ray 
albedo from the Trojan groups, which move collectively, 
or fluctuations in the SSSB 7-ray albedo at low ecliptic 
latitudes. 

In order to estimate the possible "excess" diffuse flux 
from MBAs, Trojans, and KBOs we calculated the in- 
tegrated fluxes for ecliptic latitudes |/3| < 15° and all 
ecliptic longitudes (Table [TJ "Flux" column). In order to 
increase the sensitivity, and to search for a diffuse signal 
that is centered on the ecliptic the table also includes 
fluxes for the best-fitting Gaussian centered on j3 = 0° 
and having FWHM width 12.5° ("Fitted flux" column) , 
the approximate extent of the Kuiper Belt (|Brownll2001l ). 
The fits included a constant term to account for Galac- 
tic and extragalactic diffuse emission; profiles of the fits 
are included in Figure [3 The effective PSF for EGRET 
in the 100-500 MeV range for the expected spectrum of 
the albedo emission is approximately 4° FWHM, which 
would not appreciably broaden the apparent distribution 
of 7-ray intensity. In any case, the assumption of a sin- 
gle Gaussian profile is an approximation; the contribu- 
tion from MBAs should result in an additional somewhat 
narrower but fainter component to the diffuse emission 
around the ecliptic. 

Gamma-ray emission associated with the Moon and 
Sun also contributes to the intensity of the sky at low 
ecliptic latitudes. The Moon is always within about 5° 
of the ecliptic and the profiles shown in Figure 5 un- 
doubtedly include lunar al bedo 7-ray emission, at a leve l 
of ~ 5 x 10 cm - 2 s" 1 (iMoskalenko fc Porterl l2007bD. 
As described by iMoskalenko et al.l ~i 20061 ) thesolar ra- 
diation field is a fairly bright and diffuse 7-ray source 
from inverse Compton scattering of CR electrons. The 
solar inverse Compton emission is brightest in the eclip- 
tic plane but of course depends on solar elongation angle. 
The precise contribution to the diffuse intensity at low 
ecliptic latitudes is difficult to estimate. The Sun was in 
the field of view of EGRET for only a small fraction of 
the observing time and the contribution to the total flux 
should have been less than that of the Moon. 

After the bright diffuse emission and identified point 
sources are masked from the EGRET data, no strong 
excess of diffuse emission is apparent at low ecliptic lat- 
itudes in Figure [5j The integrated fluxes are formally 
significant for the case where the Galactic plane and all 
sources identified in the 3EG catalog are masked out (Ta- 



ble [TJ Figure [5b), but the systematic uncertainties are 
comparable to the measurement. This is suggested by 
the effect on the i ntegrated flux from masking out several 
more sources that lSowards-Emmerdl (|2003l l2004h identi- 
fied as blazars (Figure [5b?). The overall average exposure 
does not change appreciably as a result of the additional 
masking but the fit flux decreased by more than 50% 
(Table [T]l. 

Based on our analysis of the EGRET data we infer 
that the diffuse emission from MBAs, Trojans, and KBOs 
has an integrated flux of less than ^6 x 10~ 6 cm" 2 s _1 
(100-500 MeV), as derived from the set of cuts 6, which 
corresponds to ~12 Lunar albedo flux units. 

5. DISCUSSION AND CONCLUSION 

The albedo 7-ray flux from MBAs can be calculated 
using eqs. (f7|)- (fTTj) and Figures [2j [3] where we assume 
that their surface material is regolith. We use the fol- 
lowing parameters: p = 2 g cm -3 for the MBA average 
density, r% — 4.565 x 10 7 cm for the radius of Ceres, 
ro = 100 cm for the smallest radius of an asteroid that 
is still an opaque target for incident CR particles. The 
central grammage in this case (ro = 100 cm) is ^400 g 
cm" 2 . Since the composition of the MBAs (and other 
SSSB populations) is mainly oxygen, this corresponds to 
^5 interaction lengths which is sufficient for the hadronic 
cascade to fully develop at the CR energies we consider. 
With these parameters, the total MBA albedo flux is 
X = Ftot/^Moon = 0.05, 0.67, 12 for extrapolation to 
small sizes with indices n = 2.5, 3.0, 3.5 (see Figure Q}, 
correspondingly. 

Similarly, for the Jovian Trojan asteroids we can esti- 
mate the 7-ray flux assuming the same size distribution 
as for MBAs, but with I ~ 5.2 AU, and r = 200 cm 
(which gives the same central grammage for p = 1 g 
cm" 3 ). We obtain X = F tot /F M oon = 0.01, 0.07, 0.8 
(averaged over their orbit) for a similar extrapolation to 
small sizes with indices n — 2.5, 3.0, 3.5. For the closest 
(4.2 AU) and the farthest (6.2 AU) distances to Earth, 
the fluxes will be 0.01, 0.1, 1.2 and 0.006, 0.05, 0.5, cor- 
respondingly. 

The KBO size distribution is known very approxi- 
mately. The second largest object of the Kuiper Belt 
after 136199 Eris is Pluto r x = 1.195 x 10 8 cm, while 
the majority of the KBOs are icy rocks and comets with 
p = 0.5 g cm -3 . To keep the same central grammage of 
the smallest body ~400 g cm -2 we have to use a larger 
minimum radius, ro = 400 cm, than for the MBAs. The 
incident spectrum of CR particles at >30 AU approaches 
the LIS spectrum which results in a factor of ^2 increase 
below ~1 GeV of the albedo flux compared to the same 
body at 1 AU (Figure [5]). For these parameters, the total 
Kuiper Belt albedo flux is X K = F^ t /F Moon = 0.2, 34, 
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1168 for n — 3.0, 3.5, 3.9, correspondingly. The removal 
of Eris and Pluto (~0.005M© combined) from the Kuipcr 
Belt and using Charon instead, r\ — 6 x 10 7 cm, results 
in the flux increase: X K = F^ t /F Moon = 0.35, 46, 1222 
for the same values of n. However, this change is simply 
the result of anchoring the power-law size distribution to 
a large body. 

Our estimates show that the albedo of MBAs and 
KBOs could account for the EGRET upper limit of the 
flux from the ecliptic. For the adopted size distributions 
of SSSBs (n = 3.0 for MBAs and n = 3.5 for KBOs), the 
KBO albedo is essentially brighter. However, if the MBA 
size distribution is somewhat steeper than our adopted 
in dex of n = 3. 0, e.g., as for the distribution proposed 
bv lChend (pool , it can account for the total albedo flux 
from the ecliptic. The SSSB 7-ray albedo, especially of 
the collectively moving Trojan groups, might be respon- 
sible for some fraction of the EGRET unidentified point 
sources at low ecliptic latitudes. 

A possible way to distinguish the albedo emission of 
MBAs and KBOs is to study the emission as a func- 
tion of solar elongation angle. In the antisolar direction, 
6 « 180°, the direction in which the Main Belt is closest 
to the Earth (~1.7 AU), the flux is predicted to be as 
much as ~5 times that in the solar direction, 6 « 0°. 
On the other hand, the brightness of the Kuiper Belt is 
only weakly dependent on the elongation angle because 
it is much further away. The positions of the Trojans on 
the sky are known, being determined relative to their re- 
spective planet (Jupiter, Neptune), making them easier 
to detect. 

The detection of the CR-induced 7-ray albedo of 
MBAs, Trojans, and the KBOs by 7-ray instruments is 
possible. At higher energies >1 GeV where the 7-ray 
albedo flux is steady and does not depend on the solar 
modulation, it can serve as a normalization point to the 
cumulative brightness of all MBAs plus KBOs. At lower 
energies <1 GeV, the component of the albedo which is 
independent of elongation, the KBO albedo, will tell us 
directly about the LIS spectrum of CRs. Therefore, the 
observations of the albedo flux can provide us with valu- 
able information about the size distributions of SSSBs 
in both regions, while the shape of the albedo spectrum 
can tell us about the LIS spectra of CR protons and he- 
lium at high energies. In turn, a detection of the MBA 
and KBO albedo at MeV-GeV energies will enable us to 
normalize properly the cumulative albedo spectrum and 
make a prediction for the intensity of the 511 keV line. 

A conservative estimate of the 511 keV line flux from 
SSSBs can be made using the upper limit derived in 
Section [U The total flux of 511 keV photons from the 
Moon is F$£ on ~ 10~ 3 Ai; rj 2.4 x 10~ 5 photons cm" 2 
s" 1 (Figure $ = MV), where AE = 0.024 MeV 
is the size of the bin containing E = 0.511 MeV. The 
total flux from the SSSBs X 511 = F^-/F^ oti can be 
calculated from eqs. (ff]l-([IT|). The SSSB albedo con- 
tribution to the 511 keV line flux within the Galac- 
tic bulge is ~ 0.72F^ 1 o 1 on , where we assumed that the 
FWHM of the bulge is ~10° (|Knodlseder et al.ll200l . 
and 20°/360° « 0.06, and we used the upper limit de- 
rived in Section |U It gives ^2 x 10 -5 photons cm~ 2 s _1 , 
which is about 2% of the total bulge emission as observed 
by the INTEGRAL (1.05±0.06) x 10" 3 photons cm~ 2 s" 1 



(jKnodlseder et apjOOg ). Since most of the INTEGRAL 
observing time was spent on observations of the Galactic 
bulge and a comparatively small fraction went into ob- 
serving regions above and below the Galactic plane, it is 
not surprising that a diffuse band near the ecliptic (the 
SSSB albedo) has not been found so far. It is interesting 
that the OSSE map of the 511 keV line has a controver- 
sial feature, the so-c alled "annihilation fountain," above 
the Galactic bulge (jPurcell et al.l Il997h which, in fact, 
may be the asteroid albedo foreground instead. Note, 
that the 7-ray spectrometer on the NEAR-Shoemaker 
spacecraf t made observatio ns of the 511 keV line from 
433 Eros (|Evans et al.l l2001) on the surface of Eros itself, 
however, it is hard to judge the absolute intensity of the 
line from the published data. 

Our estimates of the fluxes assume that the mass and 
radius distributions are valid for the whole range of 
masses, which is not necessarily true. The number of 
small bodies may be larger or smaller than the extrapo- 
lation from the distribution of more massive bodies. We 
have also assumed spherical bodies. However, the small- 
est bodies are distinctly non-spherical which would make 
them somewhat brighter than we have estimated. Thus, 
our calculations underestimate the SSSB albedo emis- 
sion. 

The bodies that are smaller than the cutoff radius 
(r = 100 cm for MBAs, 200 cm for Trojans, 400 cm for 
KBOs) will also contribute to the albedo flux. Because of 
their smaller size, only the initial stage of the CR cascade 
will develop, producing a harder albedo spectrum while 
its intensity will be reduced due to the partial conversion 
of energy of the primary CR particles into albedo 7-rays. 

We emphasise that the detection of the 7-ray albedo 
from MBAs, KBOs, and other SSSB families directly 
probes the size distribution of these bodies below the de- 
tection limit of other methods, over considerably larger 
regions of the sky. The detectability of the 7-ray emis- 
sion by these objects has implications for studies of the 
evolution of the solar and exo-solar planetary systems 
(lBrownll200l ). studies of CRs, and diffuse 7-rays. The 
GLAST Large Area Telescope (LAT) 6 , to be launched by 
NASA in May 2008, will in just one year have an essen- 
tially uniform exposure over the entire sky a factor of 40 
or more deeper than EGRET and will be free from sen- 
sitivity variations owing to ageing of consumables. This 
capability will permit detection of albedo 7-ray fluxes for 
SSSBs at even the Lunar flux level. 
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